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Abstract: In this paper, we introduce a non-minimally conformally coupled scalar field 
and dark matter in F(T) cosmology and study their dynamics. We investigate the stability 
and phase space behavior of the parameters of the scalar field by choosing an exponential 
potential and cosmologically viable form of F(T). We found that the dynamical system 
of equations admit two unstable critical points, thus no attractor solutions exist in this 
cosmology. Furthermore taking into account the scalar field mimicking as quintessence and 
phantom energy, we discuss the corresponding cosmic evolution for both small and large 
times. We investigate the cosmological implications of the model via equation of state 
and deceleration parameters of our model and show that the late time Universe will be 
dominated by phantom energy and moreover phantom crossing is possible. Our results have 
no explicit predictions for inflation and early Universe era. 
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INTRODUCTION 



Astrophysicists are convinced that the observable universe is in a phase of rapid accelerated 
expansion and commonly termed it 'dark energy' (DE) possessing negative pressure and positive 
energy density. This conclusion has been supported by several astrophysical data findings of 



supernovae SNe la 



l|, cosmic microwave background radiations via WMAP 2|, galaxy redshift 



surveys via SDSS [3] and galactic X-ray 4J. Although the phenomenon of dark energy in cosmic 
history is very recent z ~ 0.7, it has opened new areas in cosmology research. The most elegant 
and simple resolution to DE is the cosmological constant [5] but it cannot resolve fine tunning and 
cosmic coincidence problem. Hence theorists looked for other alternative models by considering 
the dynamic nature of dark energy like quintessence scalar field [(5], a phantom energy field \j\ and 
f-essence |8j]. Another interesting set of proposals to DE puzzle is the 'modified gravity' (including 
-F(T), F(R), F(G), etc) which was proposed after the failure of general relativity to explain the 
DE puzzle. This new set of gravity theories passes several solar system and astrophysical tests 

A gravitational theory can be constructed on a non-Riemannian ( Weitzenbock) manifold where 
the properties of gravity are determined through torsion of spacetime and not curvature. Some 
earlier attempts in this direction were made by Einstein himself and other researchers. A recent 
version of torsion based gravity is F(T) [l^ . where T is the torsion scalar constructed from the 



111 ] and is in good agreement with 



Numerous features of theoretical 



tetrad. Choosing F(T) = T, leads to the teleparallel gravity [1 
some standard tests of the general relativity in solar system 1(3] . 
interest have been studied in this gra v, y already include Blrkhof. theorem [13], c„ Smolog ical 
perturbations [141 ] and phantom crossing of the state parameter [15]. Moreover, the local Lorentz 
invariance is violated which henceforth leads to violation of first law of thermodynamics [ij], Q] . 
Also the entropy-area relation in this gravity takes a modified form [21[. The Hamiltonian structure 



la. The 



121- 



of F(T) gravity has been investigated and found that there are five degrees of freedom 
torsion based theory is also an alternative candidate to the mechanism of cosmic inflation 

In teleparallel gravity, the equations of motion for any geometry are exactly the same as of 
general relativity. Due to this reason, the teleparallel gravity is termed as 'teleparallel equivalent 
of general relativity'. In teleparallel gravity, the dark energy puzzle is studied by introducing a 
scalar field with a potential. If this field is minimally coupled with torsion, then this effectively 
describes quintessence dark energy. However if it is non-minimally coupled with torsion, than 
more rich dynamics of the field appears in the form of either quintessence or phantom like, or by 



experiencing a phantom crossing [18]. Xu et al 



investigated the dynamics and stability of a 
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canonical scalar field non-minimally coupled to gravity (arising from torsion). They found that 
the dynamical system has an attr actor solution and rich dynamical behavior was found. In the 
context of general relativity, a scalar field non-minimally coupled with gravity has been studied in 



311 ] . We here extend these previous studies by replacing T with an arbitrary function F(T). We 
found that such a dynamical system possess no stable equilibrium point, however rich dynamical 
behavior of quintessence and phantom energy is observed. 

We follow the following plan: In section II we write the action and equations of motion of 
our model. In section III, we give motivation for choosing particular forms of model functions. 
In section IV, we write the dynamical system in dimensionless form and discuss its stability and 
phase space behavior. In section V, we discuss the cosmological implications of our model related 
to present accelerated Universe. We provide conclusion in section VI. 

II. BASIC EQUATIONS 

We are interested in conformally invariant models i.e. models which remain invariant under 
conformal transformation. Under the rules — > Q(x)e 1 ^ (or — > Q(x)g^ u ), <f> — > Q(x)~ 1 (J), the 
equations of motion and the whole action remain invariant of a conformally invariant model. Here 
e* is the tetrad (vierbein) basis. In general for a D— dimensional gravity model, the conformal 
coupling parameter £ = 4 ^~ 2 ^ , so that in a four dimensional theory, £ = g [32|]. Further F(T) 
theories from dynamical point of view are completely different under conformal transformation 
unlike f(R) theory. Its not possible to rewrite the total action of F(T) gravity in form of teleparallel 



action plus scalar field 33J. It shows that even the pure F(T) model behaves differently under 
conformal transformation. We propose an action of F(T) gravity conformally and non-minimally 
coupled with a scalar field as 1 



d x e 



+ a<P 2 ) + \4 2 - v{cj>) + c m ] , (i) 

where e = det(e^). Here £ is a conformal coupling parameter of order unity while e = +1,-1 
represents quintessence and phantom energy respectively. The dynamical quantity of the model 
is the tetrad and scalar field eft with a scalar potential V ((/)). C m is the matter Lagrangian. It is 
assumed that both matter and scalar field are distributed as perfect fluid. The action ([I]) can be 
considered as a generalization of a teleparallel gravity non-minimally coupled with a scalar field 
[30! ] . We add some more comments on our model: the minimally coupled scalar fields are not 



1 Here 8irG = 1 
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conformal invariant, otherwise F(T) gravity is not local Lorentz invariance. After the Lorentz 
symmetry breaking, we prefer our model must have conformal symmetry. Just for pure F(T), it is 
well-known that it can be possible to write the pure F(T) action in a conformal gauge, like F(R). 
If we couple matter (here scalar field) with F(T), this matter part must have the same (conformal) 
symmetry which can be interpreted as a generalized conformal invariance. 

We assume a spatially flat Friedmann-Robertson- Walker (FRW) metric as a background space- 
time 

ds 2 = -dt 2 + a(t) 2 [dx 2 + dy 2 + dz 2 } , (2) 

where a(t) is a scale factor and = (l,a(t),a(t),a(t)). The equations of motion are obtained by 
varying the action ([I]) w.r.t. a{t) and tfi(t), we get 

6H 2 F T (1 + i4> 2 ) + ^(1 + itfW = Pm + \e4> 2 + V(4>), (3) 

4> + 3H^-e[^F-V '(0)] =0, (4) 

which are Friedmann and Klein-Gordon equations respectively. 
The second Friedmann equation is 

a 1 



8aaF T ((l + + fa 



a 4o(l + i4) 2 ){-l2F TT a 2 + a 2 Ft) 

-48(1 + ^ 2 )d A a- l F TT - a 3 ((l + ^ 2 )F + etf - 2V(^))1 . (5) 

We can rewrite ([3]) as 

3H 2 = pf + Pm , (6) 



or 



eff 

1 _ QCflf , Q QCff _ Q _ 



where p? ff is the effective energy density of scalar field written as 



+ y(^) _ _ + v ^' (2TF T - F). (8) 



Combining Eqs. © and (J7|),we get 



2« a+ H 2 = - P f, (9) 
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where p e ? is the effective pressure of scalar field given by 



„eff 



T 



8adF T ((l +^ 2 ) + £a 



6 2a(l + £0 2 )(-12F TT d 2 + o 2 F r 
-48(1 + ^ 2 )d 4 a -1 -PrT - a 3 ((l + ^ 2 )^ + e^ 2 - 2V($) 



The deceleration parameter for this model is 
a -1 



8ahF T ((l + £,(f) 2 ) + £a 



aH 2 4a(l + £0 2 )(-12F TT d 2 + a 2 F T )H 2 L 

-48(1 + C0 2 )d 4 a~ 1 F TT - o 3 ((l + ^ 2 )F + e<j> 2 - 2V(<j>)) 



(10) 



(11) 



III. CHOICE OF F(T) AND V(<j>) 



We pick a suitab 
torsion, specifically 



e f(T) expression which contains a constant, linear and a non- linear form of 
2J 



F(T) = 2<V=r + aT + C 2 , 



(12) 



where a, c and C2 are arbitrary constants 2 . The first and the third terms (excluding the middle 
term) has correspondence with the cosmological constant EoS in f(T) gravity [23]. There are 
many kinds of such these models, reconstructed from different kinds of the dark energy models. 
For example this form (I12p may be inspired from a model for dark energy from proposed form of 



the Veneziano ghost [25]. Recently Capozziello et al 26|] investigated the cosmography of F(T) 
cosmology by using data of BAO, Supernovae la and WMAP. Following their interesting results, 
we notice that if we choose 2c = \^6Ho(Q, m Q — 1), than one can estimate the parameters of this 
F(T) model as a function of Hubble parameter Hq and the cosmographic parameters and the value 
of matter density parameter. It is interesting to note that reconstruction of F(T) model according 



to holographic dark energy 



281 ] leads to the same model as (fT2|) . 



Concerning the scalar potential^we choose an exponential function which has numerous impli- 
cations in cosmological inflation 



271 ] and dark energy in the present Universe [2' 



V(<f>) = Voe^ 



(13) 



where j3 and Vq are constants. Using (|13|) . it has been shown in literature 29| that transition of 
dark energy state parameter across the cosmological boundary is possible. Also this exponential 



2 Here c does not represent the speed of light 
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FIG. 1: (Top Left) Three dimensional phase portrait of the dynamical system with e = +1. (Top Right) 
Phase space diagram of the dynamical system with e = — 1. (Bottom Left) Time evolution of dynamical 
parameters for e = +1. (Bottom Right) Time evolution of dynamical parameters for e = — 1. In the lower 
two figures, x, y, z are represented by line, dot, dot-dash respectively. 

potential is useful in assisting inflation. Without loss of generality, we assume j3 > 0. But here we 
focus only on late-time evolution of the scalar field and we do not give any prediction of our model 
for inflation. In the limit f3 — > 0, we recover the constant potential case and therefore the model is 
continuously connected with ACDM. As a remark, our choices for scalar potential and the F(T) 
can independently be responsible for DE, but in the present context of non-minimal coupling, both 
scalar field and F(T) couples non-minimally to generate the desired effect of cosmic acceleration. 
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FIG. 2: (Top Left) Time evolution of scalar fields: phantom energy (red) and quintessence (blue). (Top 
Right) Time evolution of energy densities for small time: phantom energy (blue) and quintessence (black). 
(Bottom Left) Time evolution of energy density of quintessence field for large times. (Bottom Right) Time 
evolution of energy density of phantom energy for large times. 

IV. ANALYSIS OF STABILITY IN PHASE SPACE 



We define dimensionless density parameters via 

r V 



Here x 2 and y 2 represent the density parameters of the kinetic and potential terms respectively. 
We expect interesting cases to have the scalar field rolling down the slope of the potential, as /3 > 0, 



y^. (14) 
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we should have x > 0. The equations in dimensionless variables reduce to 



dx 

In 



-2s 



L a z 



6£z 



+2£cv / 3£"V e 1/2 ^e _1 



(15) 



1 pVQy 2 

2 e 



x(-a (1 + 1/2 z 2 )-- 



-2 £ v / 6v / £~ T (a + l)xz -6ex 2 + -y 2 



dN 



+2£cx V / 3£ r Ve 1/2 V?)( a + l)-i(i + z 2 )-\ 

A ^Qy X _ y(_ a (1 + 1/2 Z 2 )- X ~- Z 2 



-n^T^ + l)xz-,ex 2 + -y 2 



+2tcx v / 3iT T yze 1/2 v?)(a + 1) _1 (1 + 



,2\-l 



6£x, 



(16) 



(17) 



where TV" = In a, is called the e-folding parameter. To discuss the stability of the system (|15D - 
([T7|), we first obtain the critical points by solving the equations (jfe = 0, 4& = 0, ^ = 0). We 
linearize the system near the critical points up to first order. After constructing a jacobian matrix 
of coefficients of linearized system, we find its eigenvalues. If all eigenvalues are negative, the 
corresponding critical point is stable (attractor), otherwise an unstable point. 



Point 


(x* , , z* ) 


Ai 


A 2 


A 3 


Stability Condition 


A 


(0,0,0) 


5+4q 


2q+1 





Unstable 


2(l+a) 


2(l+a) 


B 


(0,0,*.) 


z=(4-5£)+(5-4£) 







Unstable 


2(1 + ^)(C-1) 


2(l + ^)(?-l) 



TABLE I: Critical points and stability conditions. 



In Table I, we present the critical points, corresponding eigenvalues and stability condition. We 
notice that the system of dynamical equations admit no stable critical point. The first critical 
point A is trivial one while point B contains an undetermined component z*. For A, Ai_2 > 0, 
A3 = while for B, Ai < 0, A2 > 0, A3 = 0, thus both A and B are unstable points. 



V. COSMOLOGICAL IMPLICATIONS 



In this section, we will give some cosmological implications of our model by numerically solving 
the dynamical equations. In the top panel of figure- 1, the three dimensional phase space of x, y, z 
is plotted for two different values of e. The top left figure (for e = — 1) shows that the trajectory 
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t 



FIG. 3: Time evolution of effective EoS parameter for scalar field w c g. Solid line (dot-dash) corresponds 
to e = +1 (e = -I). 

starts from x = — 10 < causing a negative friction term —3H<f> in the Klein-Gordon equation. 
The friction term gets less dominant while scalar field and potential energy parameters dominate in 
the later evolution. It shows that the late stage evolution is determined by scalar potential alone. 
This fact is also evident from bottom left figure that only potential term dominates the dynamics 
while scalar field and the kinetic term do not contribute in cosmological dynamics. However, 
putting e = +1 in the dynamical system reverses the dynamical evolution of system as shown 
in right panel in Fig.l. The scalar potential remains vanishing while the cosmic dynamics are 
determined by kinetic term <fi and scalar field. It shows that e = +1(— 1) leads to kinetic term 
(scalar potential) dominated regimes in the late time evolution, despite the system evolves from 
same initial conditions (o(0) = 1, d(0) = H = 74.2, 0(0) = 1, 0(0) = 1, £ = 1/6). 

In Figure 2 we show the time evolution of scalar fields and logarithmic energy densities for 
e = ±1. The top left figure shows that for phantom energy, the field undergoes successive stages of 
fluctuations after a time gap of unity. Here the amplitude of phantom scalar field in each successive 
stage becomes progressively less than the previous stage. In the late time evolution, it is expected 
that the field will lose its energy and fluctuations decrease, as seen in bottom right figure. However 
for quintessence, the top left figure shows that field decays for t < 0.5, and corresponding log 
energy density stays around — 1 (p ~ e^ 1 ). 
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We define the equation of state parameter of the scalar field 




(18) 



where 




-e0 2 + V(<P) - 3H 2 [(1 + a)£<f> 2 + a], 



(19) 





-2£(1 + a)#ff - -^ 2 + - + 

-\ [pm + \4 2 + V(4>) - 3H 2 [(a + l)^ 2 + a] . 



(20) 



In figure-3, we plot the time evolution of the effective EoS parameter of the scalar field w e e for two 
available values of e. Adopting the same initial conditions as in previous figures, we observe that 
the state parameters for both models evolve from the same initial value u> e ff(0) ~ — 2 which is a 
phantom state for the state parameter. However the later evolution follows opposite trajectory i.e. 
for e = — 1, the state parameter goes on to take more negative values with time, thereby evolving 
to super-phantom state. For e = +1, the state parameter initially behaves like phantom and slowly 
evolves to cross the w e ff = — 1 boundary, while again returning to phantom regime in the late time 
evolution. 



To summarize, we investigated the stability and phase space of a non-minimally conformally 
coupled scalar field in F(T) cosmology. We found that the dynamical system of equations admit 
two unstable critical points, thus no attractor solutions exist in this cosmology. Treating the scalar 
field as quintessence and phantom energy separately, we found that for phantom energy, the late 
stage evolution is determined by the potential energy in phase space while for quintessence case, it 
is kinetic energy term that plays central role. Furthermore, the time evolution of phantom scalar 
field undergoes progressive stages of fluctuations while loosing energy density and ending up to the 
value e 4 in the relevant scale. For quintessence, the scalar field remains stable for long time with 
a constant energy density. It is interesting to note that our model correctly predicts the present 
state of the Universe is dominated by phantom energy while it will remain so in the far future. 
Further, we observed that the transition of the state parameter crosses the cosmological boundary 
twice in the case of quintessence field only. Also we noted that e = +1(— 1) leads to kinetic term 



VI. CONCLUSION 
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(scalar potential) dominated regimes in the late time evolution, despite the system evolves from 
same initial conditions. 



[1] A. G. Riess et al, Astron. J. 116, 1009 (1998); S. Perlmutter et al, Astrophys. J. 517, 565 (1999). 

[2] C. L. Bennett et al., Astrophys. J. Suppl. 148, 1 (2003). 

[3] M. Tegmark et al, Phys. Rev. D 69, 103501 (2004). 

[4] S. W. Allen et al, Mon. Not. Roy. Astron. Soc. 353, 457 (2004). 

[5] V. Sahni, A. Starobinsky, Int. J. Mod. Phy. D 9, 373 (2000); P. J. Peebles, B. Ratra, Rev. Mod. Phys. 
75, 559 (2003). 

[6] B. Ratra, P. J. E. Peebles, Phys. Rev. D 37, 3406 (1988); C. Wetterich, Nucl. Phys. B 302, 668 

(1988); A. R. Liddle, R. J. Scherrer, Phys. Rev. D 59, 023509 (1999); M.U. Farooq, M. Jamil, U. 

Debnath, Astrophys. Space Sci. 334, 243 (2011); U. Debnath, S. Chattopadhyay, I. Hussain, M. Jamil, 

R. Myrzakulov, Eur. Phys. J. C 72, 1875 (2012); K. Karami, A. Sheykhi, M. Jamil, R. Myrzakulov, S. 

Ghaffari, A. Abdolmaleki, Canadian J. Phys. 90, 473 (2012). 
[7] R. R. Caldwell, Phys. Lett. B 545, 23 (2002); R. R. Caldwell, M. Kamionkowski, N. N. Weinberg, Phys. 

Rev. Lett. 91, 071301 (2003); M. Jamil, M.A. Rashid, A. Qadir, Eur. Phys. J. C 58, 325 (2008); M. 

Jamil, Eur. Phys. J. C 62, 609 (2009); M. Jamil, A. Qadir, Gen. Rel. Grav. 43, 1069 (2011); M. Jamil, 

Int. J. Theor. Phys. 49, 144 (2010); M. Jamil, I. Hussain, M.U. Farooq, Astrophys. Space Sci. 335, 339 

(2011); M.R. Setare, M. Jamil, Gen. Relativ. Gravit. 43, 293 (2011). 
[8] M. Jamil, Y. Myrzakulov, O. Razina, R. Myrzakulov, Astrophys. Space Sci. 336, 315 (2011); M. Jamil, 

D. Momeni, N. S. Serikbayev, R. Myrzakulov, Astrophys. Space Sci. 339, 37 (2012). 
[9] S. Nojiri, S. D. Odintsov, Phys. Rept. 505, 59 (2011); S. Capozziello, M. De Laurentis, 

larXiv:1108.6266> 2; M. Jamil, D. Momeni, M. Raza, R. Myrzakulov, Eur. Phys. J. C 72, 1999 (2012); 

A. Sheykhi, K. Karami, M. Jamil, E. Kazemi, M. Haddad, Int. J. Theor. Phys. 51, 1663 (2012); S.H. 

Hendi, D. Momeni, Eur. Phys. J. C 71, 1823 (2011); A. Azadi, D. Momeni, M. Nouri-Zonoz, Phys. 

Lett. B 670, 210 (2008); D. Momeni, H. Gholizade, Int. J. Mod. Phys. D 18, 1719 (2009); M. Jamil, D. 

Momeni, F.M. Mahomed, Phys. Lett. B 702, 315 (2011); M. Jamil, D. Momeni, Chin. Phys. Lett. 28, 

099801 (2011). 

[10] K. Hayashi, T. Shirafuji, Phys. Rev. D 19, 3524 (1979); K. Hayashi, T. Shirafuji, Phys. Rev. D 24, 3312 
(1981). 

[11] F. Hehl, P. von der Heyde, G. Kerlick, Rev. Mod. Phys. 48, 393 (1976). 

[12] R. Ferraro, F. Fiorini, Phys. Rev. D 75, 084031(2007); R. Ferraro, F. Fiorini, Phys. Rev. D 78, 124019 
(2008). 

[13] X. Meng, Y. Wang, Eur. Phys. J. C 71, 1755 (2011). 
[14] R. Zheng, Q-G. Huang, JCAP 1103, 002 (2011). 



12 



[15] K. Bamba, C-Q Geng, C-C Lee, L-W Luo, JCAP 1101, 021 (2011). 
[16] M. Li, R-X Miao, Y-G Miao, JHEP 1107, 108 (2011). 
[17] N. J. Poplawski, Phys. Lett. B 694, 181 (2010). 

[18] C-Q. Geng, C-C. Lee, E. N. Saridakis, Y-P Wu, Phys. Lett. B 704, 38 (2011)4. 

[19] R. Miao, M. Li, Y. Miao, JCAP 11, 033 (2011). 

[20] B. Li, T. P. Sotiriou, J. D. Barrow, Phys. Rev. D 83, 064035 (2011). 

[21] K. Bamba, M. Jamil, D. Momeni, R. Mvrzakulov. larXiv:1202.6114l 

[22] R. M. Wald, Phys. Rev. D 48, R3427 (1993). 

[23] R. Myrzakulov, Eur. Phys. J. C 71, 1752 (2011). 

[24] M. Jamil, D. Momeni, R. Myrzakulov, Eur. Phys. J. C 72, 1959 (2012). 
[25] K. Karami, A. Abdolmaleki, larXiv:1202.2278l 

[26] S. Capozziello, V.F. Cardone, H. Farajollahi, A. Ravanpak, Phys. Rev. D 84, 043527 (2011). 

[27] A. Naruko, M. Sasaki, Prog. Theor. Phys. 121, 193 (2009); M. Bojowald, G. Calcagni, S. Tsujikawa, 

JCAP 11, 046 (2011); M. Susperregi, A. Mazumdar, Phys. Rev. D 58, 083512 (1998); A. R Liddle, A. 

Mazumdar, F. E Schunck, Phys. Rev. D 58, 061301 (1998). 
[28] M. H. Daouda, M. E. Rodrigues, M. J. S. Houndjo, Eur. Phys. J. C 72, 1893 (2012). 
[29] C. Gao, M. Kunz, A. R Liddle, D. Parkinson, Phys. Rev. D 81, 043520 (2010); B. Bozek, A. Abra- 

hamse, A. Albrecht, M. Barnard, Phys. Rev. D 77, 103504 (2008); S. Capozziello, V.F. Cardone, E. 

Piedipalumbo, C. Rubano, Class. Quant. Grav. 23, 1205 (2006); E. Elizalde, S. Nojiri, S. D. Odintsov, 

Phys. Rev. D 70, 043539 (2004); Z-K Guo, Y-S. Piao, X. Zhang, Y-Z. Zhang, Phys. Lett. B 608, 

177(2005). 

[30] C-Q Geng, C-C Lee, E. N. Saridakis, JCAP 01, 002 (2012); C. Xu, E. N. Saridakis, G. Leon, 

larXiv:1202.3781l 
[31] M. Szydlowski, O. Hrycyna, JCAP 01, 039 (2009). 
[32] G.V. Bicknell, A.H. Klotz, J. Phys. A: Math. Gen. 9, 1637 (1976). 
[33] R-J. Yang, Europhy. Lett. 93, 60001 (2011). 



